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Abstract. The network of interactions between hot gas, 
cool clouds, massive stars, and stellar remnants used in 
the chemodynamical modeling of the interstellar medium 
is investigated for the types of its solutions. In a physically 
consistent formulation for the energy transfer during con- 
densation, oscillations due to a cyclic switching between 
condensation and evaporation never occur. A closed-box 
system evolves in a hierarchy of equilibria: thermal bal- 
ance in the cloud gas, star-formation self-regulated due to 
heating of the clouds by massive stars, and the balance of 
condensation and evaporation. Except for an initial tran- 
sitory phase, the evolution of the metallicity follows that 
of the Simple Model quite closely. For galaxies with a high 
initial density, or if the condensation rate is low, the met- 
als produced by the stars may remain stored in the hot 
gas phase even in evolved systems with low gas fractions. 
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1. Introduction 

To understand the evolution of galaxies one may attempt 
to match the observational data by models which describe 
the global processes - star-formation rate, gas infall, gas 
loss - with suitable formulations. By adjusting free pa- 
rameters, quite good fits can be achieved. However, the 
number of these free parameters often is uncomfortably 
large. Moreover, this approach may not lead to a unique 
identification of the dominant physical process, as the per- 
sisting G dwarf- 'problem' (Pagel & Patchett 1975) and the 
formation of radial abundance gradients (Gotz & Koppen 
1992) illustrate. 

Our chemodynamical approach (Hensler 1987, Theis et 
al. 1992, Samland et al. 1997) tries to describe as precisely 
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as possible the known physical processes present in the in- 
terstellar medium (ISM) and its interaction with the stars. 
These local processes are coupled with the global dynam- 
ics of the gaseous and stellar components, constituting a 
physical description of the evolution of a galaxy. Since it is 
unrealistic to include all processes in their full complexity, 
one has to define a sufficiently simple but accurate net- 
work of interactions in the ISM. Our prescription, based 
on the three-component ISM model of McKee & Ostriker 
(1977) and on the formulations of Habe et al. (1981) and 
Ikcuchi et al. (1984), has successfully been coupled with 
the global dynamics for models of elliptical (Theis et al. 
1992) and disk galaxies (Samland et al. 1997). 

Another important aspect is the degree of non-linearity 
of the network which determines the behaviour of the 
model. Since the dependences of the rate coefficients are 
not well known, some caution is necessary to avoid the 
appearance of complex behaviour solely due to the mathe- 
matical formulation. We cannot yet fully settle these ques- 
tions, but what is needed is a more complete understand- 
ing of the behaviour of this type of model and an identi- 
fication of the crucial processes. 

In their chemodynamical models Theis et al. (1992) 
find that most often the star-formation rate varies slowly 
with time, but under certain conditions it undergoes 
strong nonlinear oscillations, involving the condensation 
and evaporation of the cool clouds embedded in the hot 
intercloud gas. The phases of slow variation are due to 
an equilibrium caused by the self-regulation of the star- 
formation rate (SFR) whose efficiency is reduced as the 
massive stars heat the gas by their ionizing continuum 
radiation. In a partial network with a single gas phase 
Koppen et al. (1995) show that this equilibrium results in 
a quadratic dependence of the SFR on gas density - inde- 
pendent of what was assumed for the stellar birth function. 
Under realistic conditions this is quickly reached, and it 
is unconditionally stable, quite insensitive to the rate co- 
efficients used. 
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The present study extends the network of Koppen et 
al. (1995) to two gas components, clouds and intercloud 
gas, described in Sect. 2. We investigate its behaviour by 
numerical solution which allows the extraction of analyt- 
ical conditions and relations. This permits the identifica- 
tion of the origin of the oscillations of the SFR (Sect. 3), 
and the formulation of a physically consistent description 
(Sect. 4) which leads to the identification of a second equi- 
librium, namely that between condensation and evapora- 
tion of the clouds. In Sect. 5 we extend the prescription 
to the more realistic one by Samland et al. (1997), having 
condensation and evaporation occurring simultaneously in 
a cloud population. 

2. The model for the ISM 

2.1. The description used in CDE models 

We shall consider a somewhat simplified version of the 
present CDE models which captures their characteristic 
behaviour. As in the full models, there are four compo- 
nents: the hot intercloud gas (named hereafter 'gas', with 
a mass density g), the gas in the form of clouds ('clouds' 
c), as well as massive stars (s), and low mass stars and 
remnants (r) . Between the components the following inter- 
actions are taken into account: star- formation, gas return 
from dying stars, evaporation of clouds, condensation of 
gas onto clouds, (radiative or mechanical) heating of the 
gas by massive stars, radiative cooling of the gas. The full 
network also includes other processes, such as the forma- 
tion of clouds by compression in supernova shells, dissipa- 
tion by cloud-cloud collisions. These will not be included 
in our investigation, because comparison with the results 
of the complete network showed that they do not essen- 
tially determine the type of the system's behaviour. Then 
the time evolution of the mass densities of the components 
is described by the following equations: 



9 = srj/r + Ec — Kg 

c = -\fr - Ec + Kg 

8 = £1r _ S / T 

r = (1 -^ + (l-r,)s/T 



(1) 
(2) 
(3) 
(4) 



Throughout the paper, we shall use the units parsec, 10 6 
years, and solar masses. 

Star-formation is described by the stellar birth func- 
tion used in the form of Koppen et al. (1995) 



*(c,T c ) = C n c n exp(-T c /1000K) 



(5) 



Normally we use a quadratic dependence on density (n = 
2 and Ci = 0.55). The exponential factor involving the 
temperature T c of the cloud gas describes what fraction 
of a cloud is in the form of star forming molecular clumps. 

The mass returned to the interstellar gas by dying mas- 
sive stars (with a mean life-time r — 10 Myr) is taken to 



be the fraction r/ = 0.9 of the stellar mass. Of all stars 
born, the fraction £ = 0.1 is in the form of massive stars. 

The remaining terms pertain to evaporation of clouds, 
whose rate coefficient E can be a function of densities 
and temperatures, and condensation of gas onto clouds 
(coefficient K). 

In the formulations of Hensler & Burkert (1991) and 
Theis et al. (1992) the cloudy medium is composed of 
clouds which have identical properties (radius R c , mass 
M c , density p c ) and which are embedded in the (hot) inter- 
cloud gas. One assumes pressure equilibrium (T c //i c e c = 
gT g / fig£ g ), which gives for the volume filling factor of the 
cloudy medium: 



cT c /n c 



(6) 



with the mean molecular masses (i c , fi g in the two gas 
phases. For simplicity, we set p, c — fi g — fi p equal to one 
proton mass, in all what follows. The true mean density 
p c within a single cloud is given by 



Pc = cj e c 



(7) 



Assuming that the clouds are just unstable to the Jeans 
criterion, the mass is determined by 



Mr 



~ tlr. 



knT n 



4np c \0.2Gp c 



1.5 



with Boltzmann's constant /cb and gravitational constant 
G. Hence the radius of the cloud is 



Rr. = 



kv,T r 



An 0.2Gfj, c p c 

R c [pc] = 298.0 • {TjltfK) 1 ' 2 ■ (^[Mo/pc 3 ])- 1 / 2 



(8) 



Depending on its properties, such a cloud evaporates into 
the surrounding medium due to thermal conduction from 
the ambient hot gas or the gas condenses onto it. Cowie et 
al. (1981) give a criterion for this behaviour: If the quan- 
tity defined as 



{R e \pc]y l {n B [cm- 3 ])- 1 



(9) 



1.54 10 7 K 

0.01246 ■ (T g /10 7 K) 2 • (^[pc])- 1 • ^[Mq/pc 3 ])- 1 



is smaller than 0.03, condensation occurs, otherwise the 
clouds evaporate. Note that Cowie et al. use a slightly dif- 
ferent notation: o-q. Since the clouds have identical prop- 
erties, the whole system switches between condensation 
and evaporation, depending on the criterion. In a further 
development, Samland et al. (1996) consider a cloud pop- 
ulation with a spectrum of masses, and therefore one has 
at any time clouds that evaporate as well as those con- 
densing. We discuss the implications in Sect. |H. 
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The mass loss rate per cloud due to evaporation is from 
Cowie et al.: 

|M c |[g/s] = 2.75-10 4 -T 2 - 5 - J R c [pc] (10) 
\M C \ [M /Myr] = 1.38 • 10 6 • (T g /10 7 K) 2 - 5 • i? c [pc] (11) 

For a > 1 the mass loss rate decreases with increasing a. 
In what follows, we do not consider this additional detail, 
as it alters the behaviour of the models only slightly. This 
applies also to the modifications introduced by McKee & 
Begelman (1990). From the number density of the clouds 



for convenience of analytical considerations. Likewise, the 
energy density of the cloudy medium is changed 



n c = c/M c 

one gets the rates for evaporation and condensation 

Ec = \M c \n c (13) 
Kg = \M c \n c {0m/a) (14) 

We note that with a = 0.03, on the border between evap- 
oration and condensation behaviour, the rates are equal: 



Ec = Kg 

The changes in the internal energy densities are 



(15) 



e g = h g s - g 2 A(T g ) + Ecbf c - KgbT g + bT g s/ T (16) 

with the abbreviation b = 1.5/cb/Vp = 0.013 in our units. 
The first term is the heating of the gas by massive stars. 
Because in our treatment of the stars the rate from contin- 
uous heating during the life-time of the stars hs has the 
same dependence on s as that from an explosive stellar 
death e.g. Estsis/t, processes of both types are included 
in the heating term: photoionization (with all stellar ion- 
izing photons being absorbed by the gas), deposition of 
mechanical energy by stellar winds and supernova explo- 
sions. All three processes give rate coefficients of the same 
order h ~ 200000. To permit a direct comparison with 
the model of Koppen et al. (1995) we keep h = 21900 (for 
photoionization with an efficiency of 10~ 4 for the conver- 
sion into thermal energy). A higher value gives a smaller 
self-regulated SFR and thus a longer star-formation time- 
scale. The second term is radiative cooling with a general 
cooling function A(T) (cf. Theis et al. 1992). In our units, 
and by noting the use of mass densities instead of number 
densities: 



A(T) = 1800 • vAT/IOOOOK 



(17) 



The third and fourth terms are the energy gained by the 
gas through the addition of the evaporated cloud material 
- whose temperature shall be T c - and the energy lost by 
the gas which is condensed onto the clouds. The (usually 
negligibly small) fifth term is that part of the gain if the 
stellar ejecta had gas temperature. We include this merely 



e c = h c s - c 2 A(T c ) - *&T C - EcbT c + Kgbf g 



(18) 



by gains from heating by massive stars and losses by ra- 
diative cooling, by the losses due to matter locked up into 
stars (usually a minor term), and by the losses due to 
evaporating material or by gains from incorporating the 
condensing gas. The latter temperature is designated as 
T g . For simplicity, we consider h c ~h g . 



(12) 2.2. Metallicity 



The metal masses in the gas z g = gZ g , where Z g is the 
metallicity, and the other components evolve like: 

' z 9 = s - (vz s + y^—^i 1 - z °)) + EcZ c - K 9z g (19) 



z c = -yZ c -E C Z c + KgZ g 
z s = ^Z c - sZ s /t 

Z r = (1 - 0*Z C + (1 - V)SZ S /T 



(20) 
(21) 
(22) 



The first term in the equation for z g describes the sup- 
ply to the ISM of metals from the massive stars. It is 
composed of two parts: the metals that where incorpo- 
rated into the stars at their birth and are unchanged, 
and the metals freshly synthesized in the stars. For a 
whole stellar generation the latter is ay, with the frac- 
tion a of the mass locked up into remnants which is here 
a = (1 — £) + £(1 — rf) = 1 — £77. Since the yield y (e.g. 
Koppen & Arimoto 1991) refers to the entire stellar mass 
spectrum, the massive stars alone contribute For a 
primary element such as oxygen, the yield is constant, ex- 
cept for the (negligible) factor (1 — Z s ) « 1 which takes 
into account that the primordial elements are used up with 
metal enrichment. For a secondary element (nitrogen) one 
has y oc Z s (oxygen). The abundances of every metal can 
be scaled to convenient reference values. In this paper, we 
divide the abundances by the yields, which are taken to 
be solar. 

3. Models with switching between condensation 
and evaporation 

3.1. Experiences from full numerical solutions 

This set of equations (Eqns. 1 to 22), which includes the 
switching according to the criterion of Cowie et al. (1981), 
had been incorporated into the chemodynamical models 
to compute the evolution of a closed-box system (Hensler 
& Burkcrt 1991, Theis et al. 1992), and coupled with 
1-dimensional hydrodynamics to study the evolution of 
spheroidal systems (Theis et al. 1992). It was found that 
star-formation occurs in two modes: in a continuous way 
as the consequence of self-regulation (cf. Koppen et al. 
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1995), but in a certain range of densities it would fluctu- 
ate very strongly. These fluctuations actually are regular, 
large amplitude, non-linear oscillations of the cloud tem- 
perature and hence of the star-formation rate (cf. Fig. 
7 in Theis et al. 1992). The period of the order of 100 
Myrs is longer than the cooling time scale of the cloud 
gas. This indicates that the strongly damped oscillations 
which may occur before star-formation reaches the self- 
regulated mode (Koppen et al. 1995) cannot be the ori- 
gin. Further inspection reveals that the oscillations show 
up most strongly in the cloud temperature and the gas 
density, rather weakly in the gas temperature, and were 
almost absent in the cloud density. 




-5 -4 -3-2-10 1 2 

ig(gas) 



Fig. 1. Evolution of the full system from the initial state 
(marked 'a') until the completion of the first few oscillations 
(solid line). The dashed lines descending to the right are the 
loci where the system switches from evaporation and conden- 
sation (lower curve, at point 'b') and vice versa (at point 'c'). 
The horizontal dashed line is the locus of the evaporation fun- 
nel and the dot-dashed curve depicts the condensation funnel 
(see text). 

In Fig. |l| we show the track in the (g — T c )-plane of a 
model entering the oscillatory mode and performing a few 
cycles. Starting from a point marked 'a', it quickly and 
rather directly enters a limit cycle. At point 'b' the sys- 
tem changes over from evaporation to condensation. While 
just a little amount of gas is condensed, the clouds heat 
up rapidly to a maximum, then cool along the dot-dashed 
curve, while much more gas condenses. This continues to 
point 'c' where the system changes over to evaporation. 
Cooling in the clouds follows, with initially little gas return 
which increases thereafter. Passing through a minimum 
temperature, the clouds reach a temperature marked by 
the horizontal dashed line, along which the system reaches 
again point 'b' and switches back to condensation. As 



will be shown below, the horizontal dashed line is identi- 
fied as the self-regulated star-formation mode; the passage 
through the temperature minimum is merely the transient 
while the system settles into equilibrium (Koppen et al. 
1995). 

The two dashed curves sloping down to the right delin- 
eate the border where the system changes between evap- 
oration and condensation. One notices that points 'b' and 
'c' do not lie on the same curve. This is because the gas 
temperature T g at point 'b' is about 10 percent lower than 
at point 'c'. Though this amplitude is rather modest, the 
strong dependence of the switch criterion on T g leads to a 
strong shift of the switch curve in the diagram. But these 
oscillations in T g are not essential for the behaviour of the 
model. 

3.2. Origin of the oscillations 

The fact that the numerical simulations show the appear- 
ance of the oscillations primarily in cloud temperature and 
gas density indicates that these two variables form an 'in- 
ner' system whose behaviour is governed by the values of 
the control parameters (c,T g ). A most convenient way to 
analyze this is plotting the streamlines of the equations 
in the (g — T c )-plane, with the other variables c, s, and 
T g either changing with time or held constant. We find 
that for the explanation of the fundamental structures, 
the exact dependences of the rate coefficients on density 
or temperature as given by McKee & Begelman (1990) 
or Cowie et al. (1981) are of secondary importance. It is 
completely sufficient to assume that the rate coefficients 
for condensation and evaporation are constant (K, E) . In 
the following, we regard this simplified description; for the 
criterion of switching between condensation and evapora- 
tion, the prescription to compute a is kept as given above, 
but we shall consider the threshold value of 0.03 as another 
free constant a . 

The resulting streamlines are shown in Fig. |^, which 
corresponds to Fig. [I] except for the slightly different posi- 
tion of the limit cycle. One distinguishes the two different 
regions of evaporation and condensation. In the evapora- 
tion regime, all streamlines are always oriented towards 
increasing gas density, and bunch together to form a fun- 
nel centered at constant T c w 1800 K (corresponding to 
T c = 2400 K in Fig. 0). Along this feature, the clouds 
evaporate at constant temperature, until the gas density 
exceeds a critical value - corresponding to point 'b' in the 
full model - and condensation commences. Setting K = 
in Eqns. 1, 2, and |l8| and using e c — b(cT c + cT c ), one gets 
the equation for the streamlines: 

dT c T c _ h c s - c 2 A(T c ) 

ding ® g ^ bc(sr)/T + Ec) 

At g = all streamlines are horizontal, which is merely 
a feature of the logarithmic representation and gives the 
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'g(gas) 

Fig. 2. The overall behaviour of the simplified system: stream- 
lines for the combined equations. The dot-dashed curve is the 
locus where the clouds change from condensation to evapora- 
tion, the short dashed lines mark the funnels. All streamlines 
below the horizontal dashed line point upwards. 

impression of the fan structure on the left as the lines 
bunch together into the funnel. The important horizontal 
streamline along the funnel is described by the balance be- 
tween stellar heating and radiative cooling h c s = c 2 A, i.e. 
the evolution proceeds with the clouds being essentially 
in thermal equilibrium. As these are the only features, it 
is evident that regardless of its initial state, the system 
will always end up in the funnel, along which it evolves at 
constant cloud temperature towards higher gas densities. 

In the condensation regime, all streamlines are ori- 
ented towards lower gas densities. Here, one finds a fun- 
nel, marked by the short-dashed curve sloping down from 
upper right to lower left. All streamlines, whether com- 
ing from lower or higher temperatures, merge into this 
funnel which meets the switch line at \g(g) — —2.3 and 
T c « 3000 K, and the system is changed back to evapo- 
ration corresponding to point 'c' in Fig. [IJ Setting E = 
yields the equation for the streamlines 

dT c = (h c s - c 2 A(T e )) + bKg(f g - T c ) 
ding bc(srj/T — Kg) 

Figure || depicts the whole condensation regime in the 
(g — T c )-plane. To show the behaviour at large densities, 
a lower value of T g = 20 000 K is assumed. Apart from 
a region of horizontal streamlines at g = 0, as for the 
evaporation region, one finds these structures: 

— at low gas densities, a horizontal streamline exists 
where h c s — c 2 A, the same condition as for the evapo- 
ration funnel. Its position agrees closely with the upper 
short dashed curve in Fig. 0, which is the locus where 
one has thermal balance in the clouds e c = 0. 




ig(gas) 

Fig. 3. Similar to Fig. ^, but showing the case for pure 
condensation, and assuming a temperature of the condensate 
T g = 20 000 K. On the dashed line thermal equilibrium in the 
clouds exists, the dot-dashed line shows where the cloud tem- 
perature is constant. 



— at large gas densities, a horizontal streamline appears 
where T c = f g . The dot-dashed curve (labeled 'A') 
shows where the cloud temperature does not change 
with time. The clouds are not in thermal equilibrium, 
because the condensation occurs faster than the cool- 
ing time, and the clouds are kept at the condensates' 
temperature. 

— both two parts are connected with a curve, along which 
the evolution proceeds towards thermal equilibrium in 
the clouds. Where the transition occurs, is essentially 
determined by the ratio of A and b(gK/c 2 )(T g — T c ). 

— there is an additional feature, showing up as a clump of 
streamlines: the intersection of the condensation funnel 
with the condition srj/r — Kg = 9 = at Igg = —2.6 
and T c = 1800 K (in Fig. |^) forms an attracting node. 
The evolution of the clouds proceeds in thermal equi- 
librium and in balance between stellar gas return and 
condensation. 

Independent of the initial conditions, the system always 
reaches the condensation funnel, where the clouds cool and 
the gas density decreases. However, under certain condi- 
tions, the system may enter an equilibrium evolution at 
constant gas density. 

The overall evolution can thus be pieced together: The 
two regimes of evaporation and condensation are sepa- 
rated by the curve where the criterion of Cowie et al. 
(1981) is met exactly (cr = a Q = 0.03), depicted in Fig. |[ 
One may work out how the position of this switching line 
depends on the parameters, using the prescriptions (Eqns. 
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^ to 9). For a cloud filling factor much smaller than unity, 
which is the common case, one obtains: 

a = 4.18 • 1(T 17 • T ff 2 5 T~ x g~ - 5 (23) 

The curve's location is independent of the cloud density c, 
and it is shifted towards lower gas densities, if one raises 
the threshold do or lowers the gas temperature T g , which 
is the most sensitive parameter. 

The existence of the funnels makes it relatively easy 
to predict the conditions under which oscillations occur. 
From the pattern of the streamlines it is clear that a con- 
dition for a closed track in the plane is the intersection of 
evaporation funnel and switching line should happen at 
a lower cloud temperature (and higher gas density) than 
that of condensation funnel and switch curve. The condi- 
tions for the funnels and the switch line constitute a set of 
simultaneous non- linear equations which can be solved nu- 
merically to yield the curve in the (g — T c )-plane, and the 
positions of the intersections. The stellar density can be 
estimated from the low density limit s « £t\I/ (cf. Koppen 
et al. 1995). The geometry of the funnels thus allows a 
classification of the behaviour, as a function of the 'outer' 
parameters c and T g , as shown in Fig. ^. Furthermore, the 
difference and ratio of the densities g can be used together 
with the two forms of Eqn. 1 to estimate the times spent 
in each phase, the oscillation period and the keying ra- 
tio. Because these considerations neglect the changing gas 
temperature, they predict the types of behaviour of the 
complete network only in a qualitative way. 

All these investigations - together with numerical solu- 
tions of the complete system - show that oscillations are 
restricted to a finite range in cloud density, which may 
be quite narrow (less than a decade) and which is rather 
sensitive to some of the parameters and constants in the 
model. 

During a cycle, the following sequence of physical pro- 
cesses takes place: When condensation starts, the clouds 
are rapidly heated up by the gas condensates, until the 
condensation funnel is reached. There, the clouds cool ra- 
diatively and evolve towards thermal equilibrium. When 
the switching criterion is met again, evaporation com- 
mences. Lacking the heating by the condensates, the 
clouds cool until the new thermal equilibrium is reached. 
Then the on-going evaporation of the clouds increases the 
gas density, until the condensation criterion is again ful- 
filled. 

The cycle is thus essentially driven by the deposition of 
thermal energy of the condensates into the clouds, which 
radiate the energy away in the form of line emission. Thus 
the source which makes the cycle possible is the deposition 
of sufficiently hot material onto the clouds. The tempera- 
ture T g of the condensates is an extremely sensitive param- 
eter: oscillations can only be found if the condensing gas 
remains hot T g =T g . Assumption of a lower temperature 
greatly reduces the region for the occurrence of oscilla- 
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Fig. 4. The behaviour of the simplified system as a function 
of the control parameters c and T g , as estimated with the as- 
sumption of constant stellar density: Circles denote where os- 
cillations occur (filled circles: the gas temperature is within a 
factor of 2 of the value for thermal equilibrium of the gas). 
In the white area, oscillations also occur, but with a period of 
over 5 Gyrs. Plus-signs indicate where the streamlines in both 
regimes are too fiat, and the system comes to rest at an inter- 
mediate point on the switching line. Small dots show when the 
funnels meet very closely, and the asterisks indicate that the 
system enters the 5 = equilibrium. 

tions, or can completely suppress them. This is quite eas- 
ily understood in terms of the condensation funnel (Fig. 
^|): the position of its high-density branch is determined 
by the condensate temperature. Any cooling reduces the 
range in cloud temperatures that the funnel covers and 
thus the possibility for having a cyclic path between the 
two regimes is reduced. Thus, the prescription of the trans- 
fer of energy during the condensation phase is a most crit- 
ical part of the network. 

4. The physically consistent formulation 

McKee & Cowie (1977) and McKee & Begelman (1990) 
give expressions for the mass transfer rate, but not for the 
rate of energy exchange between gas and clouds. There- 
fore one considered in the previous chemodynamical mod- 
els (e.g. Theis et al. 1992) limiting cases, such as that the 
condensing gas does not cool or cools down to some fixed 
temperature. Forming only one of the many processes of 
the ISM model network, this rather simple recipe was pre- 
ferred. Of course, it is not a physically entirely satisfactory 
description, even more so since the oscillations were found 
to be rather sensitive to the exact recipe. 

A physically consistent prescription can be formulated 
from noting that McKee & Cowie (1977) as well as Mc- 
Kee & Begelman (1990) derive their evaporation and con- 
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densation rates on the basis of a stationary solution for 
the gas flow between cloud and intercloud matter, with 
a fixed temperature profile. This means that any parcel 
of condensing gas follows the local temperature from the 
intercloud gas to the cloud interior, and that it arrives in 
the cloud with the cloud's temperature: T g — T c . As dis- 
cussed above, cooling of the condensating material results 
in a lowering of the high-density part of the condensation 
funnel (cf. Figs. | and |). If now the condensates cool 
down to the cloud temperature, the condensation funnel 
is merely a horizontal line, at the same height of the evap- 
oration funnel. Because of the opposite directions of the 
streamlines in the two regimes, the intersection of the fun- 
nels - which are characterized by thermal balance in the 
clouds - with the switching curve is an attracting node, 
in which the system will always end up. Here, the system 
evolves with the cloud gas in thermal equilibrium, and 
without any oscillations. 

The switching condition implies that Ec = Kg (Eqn. 
|l5| ) and so the evaporation and condensation terms in 
Eqns. 1, 2, [l6|, and [l8| vanish. Usually one finds that j<c, 
so that one has c « c + g, and the overall evolution is de- 
scribed by 

-^(c + g) = -f + 577/T 
s = £\J/ - s /t 
r = (1 - 0* + (1 - V)s/r 
e c = h c s - c 2 A(T c ) - *6T C 

which is nothing but the single-gas-phase system which 
has been shown to evolve almost exclusively in an equilib- 
rium between star-formation and its inhibition due to the 
heating of the ambient gas by the massive stars (Koppen 
et al. 1995). 

Note that in this more general view, the conversion 
of gaseous matter into remnants is independent of the rate 
coefficients for evaporation and condensation, and even on 
the actual prescription for the switching. But how the gas 
is divided among clouds and intercloud gas, depends on 
these specifications. 

The switching of the system between condensation and 
evaporation gives rise to non-continuous terms in the dif- 
ferential equations. However, since in both regimes the 
system evolves towards the switching condition, this ap- 
parent mathematical difficulty makes the solution instead 
more simple: It forces the system to evolve along the switch 
line. 

This permits to estimate the timescales of the equa- 
tions: For conditions appropriate to a galactic disk (c = 
0.1 M Q /pc 3 ) one has T c ss 3000 K from the figures shown. 
The gas density is typically g w 0.1. ..0.01c. The cloud fill- 
ing factor is usually quite small e c = cT c /gT g <C 1. From 
Eqn. ^3|, one gets for the temperature of the intercloud 
gas 

T g = 8.76 • 10 5 K • (ct /0.03) ' 4 • T° A g 2 w 10 7 K (24) 



which is quite insensitive to the actual gas density. The 
rate coefficients are obtained from Eqns. || to 14 

E = 1.17 ■ 10~ 17 • T 25 77 1 w 0.001 Myr" 1 

K = E--K. 0.01...0.1Myr~ 1 . 

9 

This makes these processes much slower than the cooling 
in the clouds t coo i !C = bT c /cA w 0.4 Myr, but faster than 
star-formation tsfr = rj 300 Myr. For larger heating 
coefficients than we used here, the star-formation is even 
slower (few Gyr), and the two equilibria influence each 
other even less. 

Thus, after some transient phase, the evolution pro- 
ceeds in a hierarchy of equilibria: the fast cooling in 
the clouds will establish the self-regulated mode of star- 
formation; on a longer timescale the condensation and 
evaporation will balance the distribution of the gaseous 
matter among cloud and intercloud gas; and even more 
slowly, all gas is consumed to be turned into stellar rem- 
nants. 

4-.1. Balance of condensation and evaporation 

How the gas is distributed among the gas and the cloud 
component, depends on the condensation and evapora- 
tion. As stated above, the system evolves along the switch- 
ing line, which can be formulated as: 

at Tg g T c 

Thermal equilibrium in the clouds implies constant tem- 
perature, so the last term vanishes, and with Eqn. [l^ one 

gets 

i = 5 Tg = 5 hgs - g 2 A(T g ) 
9 Tg bgTg 

The changes of the densities can be written in the low 
density limit of self-regulated star-formation s w ^rVt « 
0.03 c 2 

9 = a x 5 - - 2 (ao/0.03)-°- 4 [ a2 c 2 - a 3 5 21 (a /0.03) a2 ] 
c = _(1_ £,7)^-5 

where T g is taken from Eqn. ^4|, and using the abbrevia- 
tions 

01 = 1.79 10" 5 • (T c /3000K)-°- 4 

a 2 = 0.03 • h g « 660 

a 3 = 8.36 10 4 • (T C /3000K) 2 

One notes that the equations have a critical point which 
obeys the condition 02c 2 = a 3 g 21 (cro/0.03) - 2 , or 

9 -1/21 -2/21 
- oc c ' ct ' 
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Thus the ratio g/c due to the balance of condensation 
and evaporation increases only very slightly as the cloud 
gas is consumed. It is independent of the actual rate co- 
efficients, and only weakly dependent on the switching 
threshold. These properties are found in the numerical 
solutions, shown in Fig. After the transient evolution 
subsides - within 1 Gyr, which depends on the rate coeffi- 
cients - the system enters the equilibrium, where the above 
derived dependences are obeyed. This shows that, apart 
from transient phases, the results of the chemo-dynamical 
models are quite insensitive to the details of the recipe of 
Cowie et al. (1981), for example the correction factor <I> 
for the influence of magnetic fields. 
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Fig. 5. The mass ratio of gas and clouds, as a function of cloud 
density, for different values of the switching threshold ctq. 



The stability of the system against perturbations of the 
gas and cloud densities and temperatures are investigated 
by numerical experiments, adding or removing a substan- 
tial (50 percent) portion of the gas present in the system 
or changing strongly the temperatures. After a transitory 
phase lasting some 100 Myrs, a new equilibrium solution 
is reached which has the same character as before, i.e. 
without oscillations. In the transient there may be a few 
bumps and wiggles, but the system never breaks into os- 
cillations. This is as expected, since the disturbance moves 
the system in e.g. the (g — T' c )-plane to another position, 
from where it always evolves towards one of the funnels 
to end up in the intersection with the switch line. Since 
the character of the streamlines cannot be altered into the 
pattern responsible for oscillations, the system is also very 
robust against inflows, outflows, and external heating. 



4.2. Metallicity 

In Fig. H we depict the relation of abundance in the cloud 
component of a primary element with the gas fraction 
/gas = (c+ g) /(c+g + s + r). In the Simple Model (closed- 
box with instantaneous recycling) the gas metallicity fol- 
lows Z = —y In / gas with the true yield y, shown as the 
dashed line. For simplicity, we divide the abundances by 
the yield. The abundance rises more steeply at early times 
than the Simple Model, because we allow for a finite life- 
time of the massive stars which are responsible for the 
metal enrichment. 

For comparison, we also show the model where the 
condensing gas does not cool: the metallicity remains con- 
stant during the evaporative phases, separated by a steep 
rise when the metal-rich gas is mixed to the clouds at 
the beginning of the condensation phase. If one allows 
the gas to cool down to cloud temperature before con- 
densing, after a single evaporative intervall, the evolution 
proceeds in what seems to be the average evolution of 
the model with oscillations. It is worth emphasizing that 
for lg(— ln(/ gas )) > —0.2, the evolution of either model is 
quite close (within 0.1 dex) to that of the Simple Model, 
i.e. the effective yield is nearly equal to the true yield. 

In Fig. fj] we show the behaviour of the abundance of a 
secondary element such as nitrogen, by the N/O vs. O/H 
plot. The Simple Model predicts a strictly linear relation- 
ship, viz. a straight line with slope 1 for the logarithmic 
values. Apart from a systematically lower yield for nitro- 
gen (by 0.1 dex), the numerical models show a quite simi- 
lar slope. The oscillations are less noticeable here than in 
the abundances themselves, and in the overall evolution, 
they make the slope of the N/O-O/H relation steeper. If 
one allows for cooling of the condensing gas, the slope is 
very close to unity. 

Thus, in both aspects, the chemical enrichment of the 
clouds from which the stars are born occurs in the very 
much same way as in the Simple Model. This should not 
be too surprising, because we do consider a closed-box 
model. Moreover, the time-scales for mixing of the freshly- 
produced metals into the cloud gas are determined by 
the condensation/evaporation times, which are still short 
compared to the overall evolution, viz. the gas consump- 
tion. 



5. Models with simultaneous evaporation and con- 
densation 

The models that treat the clouds as being identical are 
certainly a simplification. For a more realistic approach, 
one would consider a cloud population with a distribu- 
tion in e.g. mass, as is done in the models of Samland et 
al. (1997). With such a description of the cloud phase, the 
Cowie et al. criterion implies that all clouds smaller than a 
certain radius will evaporate, while those larger will grow 
by condensation. Thus, one has at the same time both 
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Fig. 6. The abundance of a primary element (oxy- 
gen) in the clouds as a function of the gas fraction 
/gas = (c + g)/(c + g + s + r). Shown are models with constant 
evaporation/condensation rate coefficients, assuming that the 
condensates retain the gas temperature (stepped curve) or cool 
down to the cloud temperature. The behaviour of the Simple 
Model is shown as a dashed line. 

processes taking place, and for the overall rates one in- 
tegrates over the cloud spectrum. How this is done and 
which rates one gets, depends on any further assumptions 
on the evolution of the cloud population. Samland et al. 
obtain for the solar neighbourhood, coefficients E « 10 -3 
Myr -1 and K w 0.1 Myr -1 which remain quite constant 
throughout the evolution. 

Such a system settles quickly, within a few condensa- 
tion time-scales r con( j into the balance between condensa- 
tion and evaporation 

i = E 

c K 

which is time-independent. The overall gas consumption 
takes place much more slowly, with tsfr- For a linear SFR, 
this can easily be shown by solving the equations analyt- 
ically. Since one normally has tsfr > T cvap > T con d, most 
of the gaseous matter is in the cloud component, which is 
slowly converted into stellar remnants. 

In Fig. U we show the ratio g/c as a function of 
cloud density for models with E/K = 0.1, starting with 
an initial density c = 10 05 M Q /pc 3 . The self-regulated 
star-formation is quickly established as soon as the star- 
formation time-scale has become larger than the mean life- 
time of the massive stars (c/\& > 2t or c < 10 0,2 ). As 
long as the condensation occurs faster than the stellar gas 
return (K > 0.1 Myr -1 ), the system settles into the con- 
densation/evaporation equilibrium. In models with very 
low condensation coefficients, the ejecta from the stars 
formed in the initial period of rather high star-formation 



Fig. 7. The evolution of the abundances of a secondary (ni- 
trogen) and a primary element (oxygen), for the same models 
shown in Fig. m. 



are first accumulated in the hot gas, whose mass may re- 
main constant while the clouds continue to be consumed 
-j/cal/c. After a few condensation time-scales the gas 
condenses following the steady-state condition (Eqn. 25). 
The peak in the curve for the rather extreme model with 
K = 0.001 at lgc = — 1 marks where the gas and cloud 
phases are in equilibrium. In this model, about 10 per- 
cent of the initial mass is present in the gas component 
for about a few Gyr. 
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Fig. 8. The mass ratio of gas and clouds, as a function of cloud 
density, for different values of the condensation rate coefficient 
K (in Myr -1 ), and with the same ratio E/K — 0.1. 
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In Fig. |^ we show the ratio of the metallicities in gas 
and cloud components, to study the efficiency of the mix- 
ing between the two components. Comparing with the 
evolution of the density ratio, one notices that the metal- 
licities become equal as soon as the model reaches the 
steady-state solution (Eqn. 25). 
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Fig. 10. The oxygen metallicity as a function of the gas mass 
fraction, for the intercloud gas (and for different values of the 
evaporation rate coefficient E) and the clouds (independent of 
E). The condensation coefficient is K = 0.1, the initial cloud 
density is Co = 0.1 Mq/pc 3 . The behaviour of the Simple Model 
with a unity yield is shown as a dashed line, the short dashed 
lines for half and double yield. 



Fig. 9. Like Fig. |8|, but for the ratio of the oxygen metallicities. 
The dashed line is a model with K = 1 and E = 0.01. 

The oxygen abundance in gas and clouds as a func- 
tion of the gas fraction is depicted in Fig. [l(], compared to 
the expectation from the Simple Model. If condensation 
and evaporation occur faster than star-formation, gas and 
clouds remain well-mixed, and the model follows the Sim- 
ple Model very closely. The metallicity of the clouds is 
essentially determined by the amount of mixing due to 
condensation, it does not depend on the value of E. On 
the other hand, the gas metallicity is determined by how 
much the stellar ejecta are mixed with the low metal ma- 
terial evaporated from the clouds. For low values of E the 
mixing is poor, and the metallicity is as large as the metal 
abundance of the fresh stellar ejecta which is as large as 
2//£ ~ 10y (Eq n - 19)- Eventually, all models evolve with 
the true yield. 

In systems with a large initial density cq > 1 M /pc 3 
(shown in Fig. [ll]) the initial SFR if? w 0.03cq is greater 
than the condensation rate gK » (0.1. ..0.01 • cq) • 0.1. As 
a consequence, metals accumulate in the gas, and because 
of the small but existing reprocessing of enriched cloud 
matter the metallicity even rises. The condensation of the 
metal-rich gas occurs only at a more evolved stage, and 
then the large gas metallicity causes the clouds to have 
metallicities larger than expected from the Simple Model. 



5.1. Distribution of metal mass among the components 

In a galaxy, the total mass of a primary element such as 
oxygen is proportional to the total number of type II su- 
pernovae that have exploded until the present time, inde- 
pendent of when and where this happened. This makes the 
total metal mass an essential indicator for the global state 
of a galaxy, though it is not easy to obtain, as it requires 
the determination of not only the mean abundances of gas 
in the various forms and stars, but their mass fractions as 
well. This would involve different wavelength regions and 
abundance analysis techniques with the problems associ- 
ated combining the data. 

Is it necessary to measure all gaseous components, or 
would the observation of a single phase, e.g. the cloud 
components being observable with H II regions, suffice? 
What is the magnitude of the errors involved? In the fol- 
lowing, we use our closed-box chemodynamical models to 
provide some rough answers. Since the freshly produced 
metals are first injected into the warm or hot intercloud 
gas, observations of the dense cloud medium might miss 
an important fraction of the metals. 

In Figs. [l2| to [l4| we show how the distribution of oxy- 
gen mass over gas, clouds, and remnants depends on the 
gas mass fraction / gas , which indicates the chemical age of 
a system. The metals in chemically very young systems are 
still mainly in the hot gas. At more advanced stages, the 
metal rich gas has condensed into the clouds, which then 
contain most metal mass. In highly evolved systems, the 
metals are found in the remnants. The shape of the curves 
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Fig. 11. Same as Fig. ^Cj, but with a higher initial cloud density 
co = 10 M Q /pc 3 . 

in intermediately evolved systems, i.e. around the peak of 
the contribution from the clouds, is mainly determined by 
the ratio of the inverse time-scales for condensation and 
initial star formation K/(^>/cq) oc K ■ cq. Because in the 
closed-box model most of the star formation occurs dur- 
ing the initial phase, this time-scale determines how much 
metals are produced and which are released from the inter- 
cloud medium by condensation K. For K ■ cq < 1 M pc~ 3 
Myr" 1 (Fig. |l2|) the metals are hardly to be found in the 
gas, because of the rapid condensation into clouds. For val- 
ues smaller than about 0.1 (Fig. |l^) the metals remain in 
the gas until rather late, and because of the fast conversion 
of clouds into stars, they are swiftly locked up into rem- 
nants. As the clouds may contain less than a third of all 
the metals, measurements of the metallicity of the cloud 
gas only could be rather misleading for the determination 
of the total metal mass in a galaxy. Such a situation arises 
for low condensation rates or for high initial gas densities, 
i.e. a centrally rather concentrated proto-galactic cloud. 

The shapes of the curves are much less strongly influ- 
enced by the ratio E/K which we kept constant in the 
figures: It determines the ratio of metal mass contained 
in gas and clouds in the late evolution; thus because one 
has E < K, the clouds keep a larger portion of the metal 
mass than the gas. 

6. Conclusions 

We investigate the behaviour of one-zone chemical evo- 
lution models which utilize the chemodynamical prescrip- 
tion for a two-component interstellar medium. This is done 
by numerical solution, topological analysis of the equa- 
tions, and identification of the conditions for equilibrium 
solutions. 



Fig. 12. The fraction of the oxygen mass in the gas, clouds, 
and remnants, as a function of the gas mass fraction / ga a- 
Model parameters are Co = 0.1 Mq/pc 3 , K = 0.1 Myr -1 , and 
E/K = 0.01. 
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Fig. 13. Same as FigJT| but with K = 0.01 Myr" 1 . 

The sclf-rcgulatcd star-formation found in models with 
a single gas component (Koppen et al. 1995) is not at all 
disturbed by the additional condensation and evaporation 
processes. Due to its shorter time-scale, it dominates the 
evolution. The models follow the self-regulated solution 
very closely, without any oscillations. 

On the other hand, it is the distribution of the gas 
among the cloud and intercloud phases which is deter- 
mined by the mass exchange due to condensation and 
evaporation. The evolution also follows an equilibrium 
which is either ruled by the switching condition a = oq 
or by the relation gK = Ec. The mass of gas in the inter- 



12 



J. Koppen et al.: Condensation and evaporation of interstellar clouds in chemodynamical models of galaxies 



76% 



36% 



„ 

0) 
03 
TO 

O-0.5 

CD 
i— 



-1.5 



- 1 1 1 1 1 1 1 — 



_i i i i I i i i i_ 



-0.5 




0.5 1 

lg(-ln(f_gas)) 



Fig. 14. Same as Figjl^ but with an initial density Co = 10 
M /pc 3 . 



cloud phase is usually very small compared to that found 
in clouds. 

From realistic rate coefficients for condensation and 
evaporation in a cloud population (Samland et al. 1997) 
one finds that the evolution of the model follows a hier- 
archy of nested equilibria, which are in ascending order 
of their quite separate time-scales: thermal equilibrium in 
the cloud gas, regulation of the SFR by heating of gas due 
to radiation from massive stars, equilibrium of mass ex- 
change between gas and clouds. This is further embedded 
into the overall conversion of gas into stellar remnants. 

During the evolution in the equilibrium between con- 
densation and evaporation gK — Ec, the relation between 
metallicity and gas mass fraction follows a Simple Model 
quite closely. Gas and cloud phases are well-mixed, their 
metallicities are the same. 

Together with the quadratic dependence of the star- 
formation rate on the cloud density (cf. Koppen et al. 
1995), these aspects make up the global behaviour that is 
typical for chemo-dynamical evolution models. 

In models with low condensation coefficients and/or 
high initial cloud densities, the metal-rich gas ejected by 
the stars formed in the initial period of intense star- 
formation may remain in the intercloud medium quite 
long, before condensation into clouds commences and 
metal enrichment of the clouds occurs. 

From these models we compute how the total metal 
mass in a galaxy is distributed among gas, clouds, and 
stellar remnants and how this changes with the chemical 
age of the system. For conditions similar to the solar neigh- 
bourhood, unevolved gas-rich systems have their metals in 
both cloud and intercloud gas. In high-density models (or 
those with small condensation rates) the metals may re- 
main in the intercloud medium even in rather evolved, 



gas-rich galaxies, and those contained in the cloud phase 
may be a rather portion. Work is under way to investigate 
the properties of full 'chemodynamical' models with the 
complete network of interactions and the coupling to the 
global dynamics. 
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